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ABSTRACT 


The  distribution  of  inventory  demand  over  procurement 
leadtime  is  modeled  using  mixed  probability  distributions 
that  explicitly  account  for  the  high  incidence  of  zero 
demands  observed  in  Inventory  Control  Point  Demand  History 
Piles.  Analysis  was  limited  to  the  right-hand  tail  area  of 
the  distribution  on  the  assumption  that  that  area  is  of  crit 
ical  importance  in  determining  the  effectiveness  of  an  inven 
tory  system.  Probabilistic  models  studied  were:  1) 
Bernoulli-exponential,  2)  B  rnoulli-lognormal,  and  3) 
Bernoulli- logistic.  These  compound  distributions  were  com¬ 
pared  tc  several  standard  distributions  including  the  Pois¬ 
son,  negative  binomial  and  normal  distributions  using  a 
resampling  procedure  appropriate  in  cases  such  as  this  where 
a  paucity  of  data  exists.  Fits  obtained  from  the  75th  to 
95th  percentiles  indicated  the  mixed  distributions  may  be 
superior  as  a  class  to  the  standard  distributions  for  high- 
demand  items. 
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I.  I NTR0DPCTI3N 


Inventory  system  perfor nance1  in  a  situation  with  random 
demand  depends  upon  reorder  point2  computations.  Reorder 
point  computations  depend  upon  the  probabilistic  model  chosen 
to  represent  the  inventory  system  in  that  the  reorder  point 
is  composed  of  two  parts:  1)  The  expected  demand  during  a 
procurement  leadtime  and  2\  The  safety  level3  determined  from 
the  probabilistic  model  of  demand  over  leadtime.  Probabilis¬ 
tic  models  are  generally  only  utilized  to  represent  demand 
over  procurement  lead time  since  that  is  the  only  time  period 
when  stockouts4  potentially  occur.  The  chosen  distribution 
is  then  a  conditional  distribution  given  procurement  lead- 
time.  The  Navy  has  historically  used  three  distributions  of 
demand  for  models  utilized  at  the  Inventory  Control  Point 
(ICP)  echelon  of  the  Navy  Supply  System.  The  current  proba¬ 
bilistic  models  and  the  average  annual  demand  that  is  used  to 
determine  which  model  to  utilize  are  displayed  in  Table  1.* 


1See  Appendix  A  for  definition  of  "System  Performance". 

2See  Appendix  A  for  definition  of  "Reorder  Point". 

3See  Appendix  A  for  definition  of  "Safety  Level". 

♦See  Appendix  A  for  definition  of  "Stockouts". 

3The  Aviation  Supply  Office  (ASO)  has  recently  changed  to 
the  use  of  the  normal  distribution  as  the  model  of  choice  for 
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Table  1:  Current  Inventory  control  Point 
Probabilistic  Dei  and  Over  Leadtime  Models 

Average  Annual 

Distribution  Demand  Range 

Poisson  0-1  Low 

Negative  Binomial  1-20  Medium 
Normal  >  20  High 

If  the  actual  quarterly  demand  records  for  an  ICP  are 
examined,  one  striking  characteristic  that  is  present  in 
nearly  all  records  is  tne  high  incidence  of  zero  observa¬ 
tions.  This  is  not  surprising  when  one  considers  the  ech¬ 
eloning  of  the  supply  system  and  the  fact  that  users  only 
place  a  demand  on  the  wholesale  system  as  their  reorder 
points  are  reached.  Since  most  activities  stock  at  least  one 
quarter’s  expected  demand  for  an  item,  reorders  are  expected 
at  most  four  times  a  year.  By  the  natural  phase  differences 
of  each  activity's  reorder  actions,  the  ICP  may  experience 
zero  demands  during  a  given  quarter  for  even  high-demand 
items,  and  will  surely  often  experience  quarters  of  zero 
demands  for  medium-  and  low-demand  items.  A  previous  study 
of  Pacific  Fleet  Combat  Stores  Ships'  demand  data  [Ref.  1] 
found  that  the  frequency  of  zero  observations  in  those  demand 


all  demand  categories.  Ships  Parts  Control  Center  (SPCC) 
still  uses  all  the  distributions  in  Table  1. 


records  appeared  to  have  a  low  correlation  with  the  overall 


demand  level  (i.e.  high,  medium  or  low).* 

A  desirable  model  for  an  inventory  system  should  compute 
reorder  levels  that  accurately  correspond  to  a  specified 
stockout  risk  based  on  the  level  of  actual  or  anticipated 
demand  for  each  item.  If  such  a  model  were  available,  the 
inventory  system  could  avoid  the  added  expense  of  either  hav¬ 
ing  too  much  or  too  little  stock.  One  problem  with  the  dis¬ 
tributions  in  Table  1  is  that  none  can  account  for  the 
probability  mass  at  zero;  therefore,  none  can  accurately  com¬ 
pute  reorder  levels. 

This  study  explores  the  use  of  compound  probability  mod¬ 
els  and  their  potential  to  more  accurately  achieve  the  goals 
presented  above  by  explicitly  accounting  for  the  probability 
mass  at  zero.  The  results  obtained  should  be  considered  as  a 
first  step  in  the  exploration  of  a  class  of  models  which  have 
heretofore  received  little  attention  in  the  context  of  model¬ 
ing  inventory  demand. 


‘See  Section  II.B.1  for  a  more  detailed  discussion  of 
this  study. 


II .  BACKGROUND 


A.  CURRENT  MODEL  SELECTION  TECHNIQUES 

In  the  past,  the  usual  approach  of  choosing  the  probabi¬ 
listic  model  for  lead  time  demand  was  based  on  matching  the 
empirical  cumulative  distribution  function  determined  from 
the  actual  demand  to  various  theoretical  cumulative  distribu¬ 
tion  functions.  The  shapes  of  the  two  curves  are  compared 
using  the  Kolmogorov-Smirno v  Statistical  Goodness  of  Pit 
test.7  The  Chi-Squared  test  has  historically  not  been  used 
because  it  is  not  very  powerful  when  the  number  of  observa¬ 
tions  is  small  (supply  demand  data  is  retained  for  at  most 
twelve  quarters  and  usually  only  eight  guarters  at  the  ICP 
level).  The  Kolmogorov-Smi raov  test  provides  only  a  relative 
measure  o£  goodness  of  fit,  but  it  is  usable  with  very  small 
sample  sizes.  The  goal  then  is  to  select  the  model  with  the 
best  agreement  between  the  empirical  and  theoretical  cumula¬ 
tive  distribution  function  curves. 

The  Navy  Fleet  Material  Support  Office  (FMSO)  has  uti¬ 
lized  the  above  Kolmogorov- Smirnov  test  to  evaluate  several 
theoretical  distributions.  Host  recently  [Ref.  2]  FHSO 


7 A  good  description  of  the  Kolmogorov-Smirnov  one  and 
two-sample  tests  may  be  found  in  Siegel,  s.,  Mon  Parametric 
Statistics,  p.47-52  and  127-136,  HcGraw  Hill,  1956. 
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evaluated  seven  distributions8  at  the  90*  confidence  level 
using  demand  data  from  NAS  Brunswick.  The  use  of  demand  data 
from  a  stock  point  vice  an  ICP  assumes  the  underlying  demand 
behavior  at  the  stock  point  and  ICP  are  similar.  Though  this 
was  never  validated  directly,  the  NAS  Brunswick  demand  data 
were  deemed  to  have  " unique  qualities  better  suited  for 
demand  analysis"  as  follows: 

1.  Twelve  quarters  of  demand  and  demand  frequency  data 
were  available  per  item,  50*  more  than  available  in  ICP 
files, 

2.  New  item  identifiers  that  permit  the  selection  of  only 
steady  state  demand  items  were  available, 

3.  Reliability  of  the  data  was  established  from  prior  FHSO 
studies. 

The  FMSO  analysis  showed  that  Navy  demand  patterns  are 
very  poorly  modeled  by  any  of  the  seven  distributions  stud¬ 
ied.  No  new  models  were  proposed  and  no  attributes  of  the 
demand  patterns  were  indicated  as  the  most  significant  factor 
in  the  failure  of  the  standard  distributions  to  model  Navy 
demand  patterns. 

There  are  several  problems  with  the  way  the  Kolmogorov- 
Smirncv  test  was  used  in  the  FNSO  study: 

•Distributions  tested  were:  Normal,  Negative-Binomial, 
Poisson,  Logistic,  LaPlace,  Gamma  and  Uniform. 


13 


1.  The  procedure  is  not  strictly  applicable  in  the  case  of 
discrete  distributions  such  as  the  Poisson  and  negative 
binomial. 

2.  The  test  assumes  that  the  parameters  of  the  distribu¬ 
tion  being  tested  are  completely  specified  in  advance, 
i.e.  not  estimated  from  the  data.  (Lilliefors  [Ref.  3] 
developed  alternate  tables  for  use  in  testing  the  expo¬ 
nential  distribution  when  parameters  are  estimated  from 
the  data.  This  procedure  was  used  in  Reference  1.) 

3.  The  Kolmogorov- Smirnov  test  evaluates  goodness-of-fit 
over  the  entire  range  of  the  distribution.  For  inven¬ 
tory  problems,  the  region  of  greatest  interest  is  the 
right-hand  tail  of  the  distribution. 

B.  PREVIOUS  RESEARCH 

Research  into  distributions  which  have  properties  that 
more  closely  model  actual  demand  patterns  has  been  very  lim¬ 
ited  and,  when  conducted,  often  did  not  explore  anything  but 
the  usual  standard  distributions.  Two  projects  though  stand 
out  because  of  their  innovative  use  of  non-standard  distribu¬ 
tions  and  the  good  fit  achieved. 

1  •  Beopulli-esB.onenti, al  Distribution 

In  Reference  1,  Pacific  Fleet  Combat  stores  ship 
demand  was  modeled  using  the  mixed  Bernoulli-exponential 


distribution.  This  distribution  was  chosen  because  it  was 
observed  that  nearly  all  line  items  had  an  unusually  high 
occurence  of  demands  of  zero.  Analysis  of  the  data  indicated 
that  with  high  probability,  the  number  of  non-zero  observa¬ 
tions  for  a  given  item  was  unrelated  to  the  average  value  of 
those  non-zero  observations.  The  conclusion  was  that  the 
demand  process  could  be  thought  of  as  two  independent  subproc 
esses  with  one  process  determining  whether  a  demand  would 
occur  or  not  and  the  other  determining  the  quantity  of  the 
demand,  given  the  demand  did  occur.  The  former  process  was 
modeled  as  a  Bernoulli  process  with  parameter  p,  p  being  the 
probability  that  a  demand  did  occur.  The  latter  process  was 
modeled  as  a  continuous  exponential  process  based  on  explora¬ 
tory  data  analysis.  The  resulting  probability  distribution 
function  for  the  number  of  units  demanded  in  a  leadtime  for 
the  Bernoulli -exponential  distribution  is  given  by: 

1  -  p  ,  x  =  0  (1) 

h(x)  =  < 

pAex  p  (- Ax) ,  x  >  0, 

\ 

where: 

1/A  *  the  expected  value  of  demand,  given  that  the 
demand  is  greater  than  zero, 
p  *  probability  that  demand  does  occur. 


The  complementary  cumulative  distribution  function  is  given 

t 

by: 


H  (x)  *  p  exp  (-Ax)  ,  x  >  0,  (2) 

where  H  (x)  is  the  probability  that  demand  will  not  exceed  x; 
this  is  equivalent  to  the  complementary  cumulative  distribu¬ 
tion  function  used  by  Hadley  and  Whitin  [Ref.  4]. 

The  hypothesis  that  the  demand  of  Combat  Stores  Ships 
could  be  modeled  according  to  the  above  probability  distribu¬ 
tion  was  tested  using  the  X olmogorov- Smirnov  test  described 
earlier  plus  tests  on  the  theoretical  risk  and  the  observed 
risk  for  values  of  risk  in  the  upper  right  hand  tail  area  of 
the  distribution.  The  results  of  these  tests  on  five  inde¬ 
pendent  samples  from  the  demand  data  available  provided 
strong  evidence  that  the  Bernoulli-exponential  distribution 
describes  the  demand  for  any  given  stock  item  very  well. 

2 •  Loaarithmic-P oisso a -Gamma  Distribution 

K  model  for  the  distribution  of  demand  during  pro¬ 
curement  leadtime  [Ref.  5]  was  derived  under  the  following 
assumptions: 

a.  Requisitions  occur  according  to  a  stationary  Poisson 
Process, 

b.  Requisition  sizes  follow  a  logarithmic  distribution. 
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c.  Leadtine  is  a  random  variable  with  the  gamma 


distribution. 

The  resulting  probability  function  for  the  number  of  units 
demanded  in  a  leadtime  is: 


x  ,  k 

.k 


h(x) 


k=l 


and 
h(0)  = 
where: 


3=1 

for  x  *  1,  2#  3 ,  ... 


f  8  ' 


a, 0  =  parameters  of  the  Gamma  distribution  defining 
leadtime, 

q  *  parameter  of  the  logarithmic  distribution 
solved  for  by  interval  bisection  from 
E(X)  «  -0  /  (1  -  0 )  In  {1  -  0)  , 

X  *  requisition  arrival  rate, 
c  =  -  A  /  In  (1  -  0) 

Sx  k =  Stirling  numbers  of  the  first  kind  computed 


from  the  recursion 

sx,k *  sx-i, k- 1  +  “  1)  ^x-l,k 

for  k  =  1,  2,  ...»  x  and  x  =  1,  2,  ... 

with  Sv  n  -  3  for  all  x 


This  distribution  is  called  the  Logarithmic-Poisson-Gamma 
(LPG)  .  It  is  defined  froa  four  paraaeters  a ,  8,  0,  and  X 
(0  and  X  deteraine  c)  . 

For  more  efficient  computations,  a  recursion  equation 

for  integer  a  was  developed. 

x 


h(x) 


(4) 


k=l 


where: 


0 

Tx,k  =  “  Cc2(ct  *  k  ~  1>rx-l,k-l  +  (*  - 
C2  =  c  /  (X  ♦  8) 

and  0,  c,  a,  8  and  X  defined  as  in  (3). 

Computing  reorder  points  from  fractiles  of  the  LPG 
distribution,  even  using  the  simplified  version  of  (4) ,  is 
generally  too  time  consuming  for  many  real  applications,  as 
a  result,  an  approximation  was  developed  and  tested  that  used 
a  scaled  version  of  the  Poisson  distribution  to  approximate 
the  negative  binomial  distribution.  The  resulting  approxima¬ 
tion  is: 


P  (Z  (t)  *  lex}  * 


(a  ♦  x  -  1)  ! 
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(5) 


x!  (  a  -  1) 


\v  +  8/\y  +  6/ 

for  x  *  0,  1,  2,  ... 


where : 


Z(t) 


number  of  units  demanded  in  time  t 


a,  8,  c  and  0  defined  as  in  (3). 

The  authors  conducted  subjective  tests  of  goodness  of 
fit  to  their  LPG  model  and  its  derivatives  using  Air  Force 
consumable  items  and  concluded,  using  as  an  example  one  item, 
that  there  was  a  “very  close  agreement  between  the  observed 
and  the  predicted  cumulative  distribution  functions  for  this 
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III.  PROPOSED  MODELS 


A.  COMPOUND  BERNOULLI  MODELS  IN  GENERAL 

The  results  of  the  previous  research  indicate  that  a 
potential  exists  for  improving  inventory  performance  by  the 
use  of  slightly  more  complex  models.  The  focus  of  this  study 
is  on  the  compound  Bernoulli  models,  leaving  further  analysis 
of  the  LPG  model  and  its  derivatives  to  others.  The  use  of 
compound  Bernoulli  models  has  several  advantages  over  stan¬ 
dard  simple  distributions: 

1.  The  model,  by  necessity,  has  an  added  parameter  whose 
sole  function  is  to  estimate  the  probability  of  a  non¬ 
zero  demand. 

2.  For  simple  models  defined  only  over  the  positive  real 
numbers,  the  addition  of  the  Bernoulli  parameter  simply 
results  in  the  scaling  of  the  cumulative  distribution, 
thus  adding  little  additional  complexity  to  the  model. 

3.  By  explicitly  accounting  for  zero  demands,  compound 
models  account  for  more  of  the  observed  variance  than 
do  simple  models. 

Compound  models  are  not  without  disadvantages: 

1.  Each  parameter  added  to  the  model  must  be  estimated  for 
each  individual  stock  item.  For  a  large  inventory 


system  such  as  the  Navy's,  this  equates  to  considerabl 
additional  computational  requirements. 

2.  Additional  parameters  must  be  stored  for  future 
retrieval  or  else  computed  from  stored  data  when 
required,  thus  either  extending  computational  time  or 
else  requiring  additional  online  storage  space. 

The  Bernoulli  compound  models  all  have  similarly  shaped 
cumulative  distribution  functions.  Pigure  1  is  a  sample  dis 
tribution  function  for  the  compound  Bernoulli-lognormal  dis¬ 
tribution  with  a  Bernoulli  factor  of  p  equal  to  0.6,  a  mean 
of  1.0  and  a  standard  deviation  of  0.5.  Note  that  the  dis¬ 
tribution  has  a  mass  of  zero  for  negative  observations,  a 
mass  of  1-p  at  zero  and  the  usual  cumulative  distribution 
function  shape  for  positive  observations  but  of  mass  p  vice 
mass  one. 

B.  HODEL  DESCRIPTIONS 

Three  compound  Bernoulli  models  were  formulated  for  eval 
uation.  Each  model  was  derived  from  the  base  distribution 
after  compounding  with  the  Bernoulli  process.  An  additional 
requirement  was  made  that  the  inverse  of  the  compound  model 
must  be  computationally  relatively  simple. 

1.  ai  lice  i£  cue  nti^l_D  1st  ribut  ion 

This  distribution  was  derived  in  Reference  1  and  the 
derivation  will  not  be  repeated  here.  Equations  (1)  and  (2) 
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Figure  1:  Compound  Cumulative  Distribution  Function, 
Bernoulli-lognormal  Distribution 

are  the  probability  density  and  complementary  cumulative 
distribution  functions  respectively.  The  Bernoulli- 
exponential  inverse  complementary  cumulative  distribution 
function  is  computed  as: 

/ 

0  ,  H(x)  >  p 

x  =  * 

M[ln(p)  -  In  (H  (x ) )  ],  H(x)  <  p 

\ 

The  parameters  of  this  distribution  are  estimated  from  the 
sample  data  as  follows: 


p  =  number  of  non-zero  demand  observations  in  sample 
divided  by  sample  size, 
y  *  x  /  p,  where  x  is  the  entire  sample  mean, 
or 

y  =  x,  where  i  is  the  sample  mean  of  only  the  non-zero 
demand  observations. 

2.  Bernoulli-lognormal,  Distribution 

The  lognormal  distribution  was  chosen  as  a  candidate 
for  testing  because  of  several  desirable  properties: 

a.  It  is  defined  for  all  real  numbers  greater  than  zero, 

b.  Its  inverse  is  readily  computed  using  inverse  normal 
approximations  that  are  well  documented  for  accuracy, 

c.  Alone,  the  lognormal  distribution  does  not  allow  the 
case  of  zero  demand,  but  when  combined  with  the  Ber¬ 
noulli  distribution,  the  compound  Bernoulli-lognormal 
distribution  defines  all  the  expected  demand  values 
greater  than  cr  equal  to  zero. 

utilizing  the  Bernoulli  parameter  p  as  a  scaling  fac 


ter  for  the  lognormal  density  function,  the  compound 
Bernoulli- log normal  density  function  may  be  expressed  as: 


p  =  probability  of  a  non-zero  demand, 

U  =  mean  of  the  natural  log  of  all  demands  greater  than 
zero, 

o  =  standard  deviation  of  the  natural  log  of  all  demands 
greater  than  zero, 

<p  =  Standard  normal  density  function. 

Integrating  (7)  results  in  the  complementary  cumulative  dis¬ 
tribution  function  given  by: 

P  ,  X  =  0  (8) 

H  (X)  = 

r  /in(x)  -  uv 

p[,  '  ♦( - 5"”  J  1  >  0 

where : 

$  *  Standard  normal  cumulative  distribution  function. 

The  Bernoulli-lognormal  inverse  cumulative  distribution  func¬ 
tion  is  easily  derived  from  (8)  as: 

0  ,  H(x)  >  p  (9) 

exp[y  ♦  a$-i  (1  -  H(x)/p)  ],  H  (x)  <  p 

The  parameters  are  estimated  from  the  sample  data  in  the 
usual  manner:  . 

p  =  number  of  non-zero  demand  observations  in  sample 
divided  by  sample  size. 


y  =  x,  where  x  is  the  sample  mean  of  the  natural  log  of 
all  non-zero  demand  observations, 

A 

o  -  s,  where  s  is  the  sample  standard  deviation  of  the 
natural  log  of  all  non-zero  demand  observations.  An 
alternate  method  for  estimating  the  mean  and  stan¬ 
dard  deviation  of  the  lognormal  distribution  is  to 
use  the  method  of  moments.  The  procedure  in  effect, 
transforms  the  mean  and  standard  deviation  obtained 
from  the  untransformed  sample  data  vice  transforming 
the  data  first  and  then  computing  the  sample  mean 
and  standard  deviation.  This  method  involves  fewer 
logarithms  and  will  run  faster  on  a  computer.  The 
procedure  is  described  in  Appendix  D. 

3.  Bernoulli-logistic  Distribution 

The  logistic  distribution  is  a  pseudo-normal  distri¬ 
bution  that  is  similar  in  shape  to  the  normal,  but  is  easier 
to  handle  mathematically.  Its  density  function  is  defined  as 
[Ref.  6]s 


r  / tt(x  -  u  )\  (10) 

f(x)  *  — --  sech *( - -r — ),  -»  <  x  <  », 

4 \J%o  \  2  JTo  J 


while  its  cumulative  distribution  function  is  defined  as: 


i  r  /"-<x  -  m  ) \ i 

'(x)  =  -  1  ♦  tanhl - 1  ,  - 00  <  x  < 


where: 


V  =  mean  of  the  entire  sample, 
o  *  standard  deviation  of  the  entire  sample. 

With  the  logistic  distribution,  it  is  no  longer  a  simple  mat¬ 
ter  of  scaling  the  distribution  by  the  Bernoulli  p  parameter, 
since  this  distribution  is  defined  over  the  entire  real  line 
and  not  just  on  the  positive  half.  Consequently,  the  distri¬ 
bution  requires  a  new  constant  of  integration  to  replace  the 
1/2  used  in  (11).  Integrating  (10)  from  0  to  *  provides  the 
new  constant  as: 


♦  tanh 


(12) 


which,  when  compounded  with  the  Bernoulli  distribution,  gives 
the  Bernoulli-logistic  distribution  density  function  as: 


/  1  -  P 


h  (x)  = 


,  x  =  0  (13) 


1  ♦  tanhj - -  ). 

Wn/j 


/7T(X  -  U  )\ 

sech2( - )  ,  x  >  0 

\  2\[3o  J 


,  otherwise? 


where : 


u  =  mean  of  all  non-zero  demand  observations. 


cr  =  standard  deviation  of  all  non-zero  demand 
observations . 

After  integrating  (13),  the  complementary  cumulative  distri¬ 


bution  function  may  be  , writ  ten  as: 


H(x)  *  P<1  -  - - - - - )rtanh 


1  ♦  tanh 
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‘C'J'o) 


♦  tanh 


ff(x  -  uM  (14) 


2  v/To 


for  x  >  0 


The  inverse  cumulative  distribution  is: 


H(x)  >  p  (15) 


I - 

u  ♦ - In/  r 

T  H  (x) 


,  *  ” r"m-<  p 


2)/?  a 


The  parameters  are  estimated  from  tha  sample  data  in  the 
usual  manner: 

p  *  number  of  non-zeco  demand  observations  in  sample 
divided  by  sample  size, 

y  *  x,  where  x  is  the  sample  mean  of  all  non-zero  demand 
observations , 

a  =  s,  where  s  is  the  sample  standard  deviation  of  all 
non-zero  demand  observations. 
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17.  EVALPATION  PBOCEDP8E 


A.  THE  DATA 

In  order  to  test  the  validity  of  the  proposed  models, 
samples  of  actual  demand  were  obtained  from  the  Operations 
Analysis  Department  at  FMSO.  The  data,  accumulated  from  the 
demand  history  files  of  the  Aviation  Supply  Office,  was  orig 
inally  used  as  input  data  for  the  5A  (Aviation  Afloat  and 
Ashore  Allowance  Analyzer)  [Bef.  7].  The  data  consists  of 
1587  consumable  IB  cog  items  and  2892  non-program-related* 
repairable  2B  cog  items.  The  information  for  -each  item  is 
contained  in  one  master  record  and  several  subrecords.  The 
master  record  contains  identifying  information  on  each  item 
such  as  the  national  stock  number,  replacement  price,  etc. 
Each  subrecord  contains  up  to  forty-six  demand  records,  each 
including  the  demand  guantity  and  the  day  of  the  demand.  A 
complete  record  layout  is  contained  in  Appendix  E.  To  ease 
processing,  the  Julian  Dates  of  the  original  demands  were 
replaced  by  FHSO  with  a  sequential  date  ranging  from  1  to 
1500  representing  approximately  four  years  of  available 
history. 


♦See  Appendix  A  for  definition  of  "program  related" 
items. 
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For  this  application,  additional  screening  was  conducted 
to  group  the  desands  into  thirty-day  buckets,  thus  providing 
a  desand  tiae  series  of  forty-eight  observations.  Further 
editing  vas  inquired  prior  to  analysis  to  reaove  the  negative 
desands  noted  for  sany  of  the  itess.  It  was  assused  that  the 
negative  observations  were  the  result  of  cancellations  of 
previous  desands  which  were  never  filled.  Therefore,  when¬ 
ever  a  negative  desand  was  encountered,  the  desand  series  was 
searched  backwards  for  the  first  positive  desand  equal  to  or 
larger  than  the  negative  desand.  If  such  a  desand  was  found, 
the  positive  desand  quantity  was  reduced  by  the  absolute 
value  of  the  negative  quantity. and  the  negative  quantity  was 
set  equal  to  zero.  If  no  offsetting  positive  desand  was 
found,  the  negative  desand  quantity  was  still  set  equal  to 
zero.  The  resulting  edited  desand  tine  series  was  used  as 
the  base  for  all  further  analysis.  Two  additional  demand 
series  were  created  fron  the  edited  series: 

1.  A  desand  series  of  only  the  non-zero  demands, 

2.  A  demand  series  of  the  log  of  the  non-zero  demands. 

The  sample  sean  and  sanple  standard  deviation  for  each  of  the 
three  demand  series  were  computed10  and  the  group,  represent¬ 
ing  one  line  item,  was  accepted  for  futher  processing  if: 


*■ °Stationarity  was  assumed  for  the  demand  time  series 
under  study.  If  a  trend  were  in  fact  present,  the  resampling 
technique  effectively  would  eliminate  it  by  shuffling  the 
observed  demands  into  random  order. 


1.  There  were  at  least  two  non-zero  demands,  and 

2.  The  standard  deviation  of  all  three  series  was  non¬ 
zero. 

Once  accepted,  the  sample  mean  and/or  sample  standard 
deviation  of  the  appropriate  series  was  used  to  estimate  the 
value  of  x,  such  that  for  a  given  probability,  demand  will 
not  exceed  that  value  of  x.  The  values  of  x  were  computed 
for  each  model  tested  at  several  probabilities  from  the 
inverse  cumulative  probability  functions.  Table  2  displays 
the  sample  mean  and/or  sample  standard  deviation  required  for 
each  model  investigated. 

B.  THE  RESAMPLING  PROCEDOR E 

The  resampling  procedure  used  was  chosen  over  the  more 
traditional  methods  discussed  in  Chapter  II  because  it  pro¬ 
vides  a  method  of  comparing  a  theoretical  distribution  to  the 
sample  data  at  specified  probabilities  and  is  applicable  in 
cases  such  as  this  where  there  is  a  paucity  of  data  available 
for  analysis.  This  procedure  is  similar  to  the  "bootstrap” 
procedure  used  by  Efron  [Ref.  8].  The  idea  behind  this  pro¬ 
cedure  is  tc  randomly  sample  with  replacement  the  series  of 
available  data  to  create  additional  pseudo-samples  that  pos¬ 
sess  the  same  statistical  properties  as  the  original  sample. 
The  desired  statistical  property  can  then  be  estimated  from 


Table  2:  Required  Data  for  Inverse  computations 


MODEL  MEAN 

1.  Exponential  A 

2.  Normal  A 

3.  Poisson  A 

4.  Negative  Binomial  A 

5.  Lognormal  C 

6.  Logistic  A 

7.  La  Place  A 

8.  Bernoulli-exponential  B 

9.  Bernoulli- lognormal  C 

10.  Bernoulli-logistic  B 

w  her  e: 


STANDARD 

DEVIATION 

N/R 

A 

N/R 

A 

C 

A 

A 

N/R 

C 

B 


A  is  from  edited  demand  series 

B  is  from  non-zero  demand  series 

C  is  from  log  of  aon-zero  demand  series 

N/R  indicates  parameter  Not  Required  for  this  model 


each  pseudo  sample.  With  a  sufficient  number  of  repetitions, 
the  distribution  of  the  estimated  property  is  known  to  be 
normal  from  the  central  limit  theorem  with  mean  equal  to  the 
population  mean  and  standard  deviation  equal  to  the  popula¬ 
tion  standard  deviation  divided  by  the  number  of  repetitions. 

Thus  by  the  use  of  the  resampling  procedure,  it  is  possi¬ 
ble  to  study  theoretical  distributions  at  points  in  their 
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right  hand  tail  area,  which  is  the  region  of  most  interest  in 
computing  safety  levels  in  inventory  models.  The  resampling 
procedure,  as  applied  in  this  study,  is  a  relatively  straight 
forward  application  of  the  following  technique: 

1.  The  edited  series  for  each  sample  of  n  demand  observa¬ 
tions  is  treated  as  the  population  from  which  the  ran¬ 
dom  samples  are  drawn. 

2.  A  random  sample  of  size  n  is  drawn  from  this  popula¬ 
tion,  creating  a  pseudo  sample.  The  pseudo  sample  is 
not  a  permutation  of  the  original  population  since  the 
sample  values  are  selected  with  replacement  from  the 
original  population.  This  is  easily  accomplished  in  a 
computer  by  generating  n  uniform  random  integers  in  the 
range  from  1  to  n,  and  using  these  numbers  as  sub¬ 
scripts  to  select  the  sample  from  the  original 
population. 

3.  The  pseudo  sample  percentile  is  computed  for  each  rep¬ 
etition  from  (16)  and  is  compared  to  the  theoretical 
percentile. 

p  *  (number  of  demand  observations  <  x)  /  n  (16) 

The  various  values  of  x  are  computed  from  the  inverse 
cumulative  probability  functions  for  each  model 


evaluated 


Under  the  null  hypothesis  that  the  sample  demands  come 
from  a  particular  probability  distribution,  the  expected 
value  of  p  should  equal  p,  the  theoretical  percentile: 

E(  p]  =  p  under  HQ  (17) 

and  the  distribution  of  p  should  be  Normal  (p,a2/n) .  The  per¬ 
centile  p  is  estimated  as  above  and  a2  is  estimated  from 
standard  binomial  results  as: 

Var[p]  *  p(1  -  p)  /  a.  (18) 


C.  MEASURES  OF  EFFECTIVENESS 

When  the  above  precedure  is  executed  at  several  theoreti¬ 
cal  percentiles,  the  fit  of  the  sample  data  may  be  compared 
at  several  locations  in  the  theoretical  distribution  and  the 
overall  fit  evaluated.  Several  common  measures  of  effective¬ 
ness  are  available  to  provide  a  criteria  on  which  to  judge 
the  success  or  failure  of  a  model.  Some  of  those  available 
are : 

1.  Maximum  absolute  deviation, 

2.  Mean  squared  error, 

3.  Algebraic  sum  of  errors. 

Method  three  was  eliminated  since  errors  may  offset  one 
another  resulting  in  a  seemingly  good  fit,  but  in  reality  a 


very  poor  fit.  Method  one  is  potentially  a  good  way  to  eval¬ 
uate  error,  but  is  generally  not  as  sensitive  as  method  two. 
Therefore,  the  aean  squared  error  of  the  pseudo  saaple  per¬ 
centile  estimates  was  chosen  as  the  measure  of  effectiveness 
of -'choice  and  was  accuaulated  as: 


i=l 


D.  COMPOTES  VALIDATION 

The  computer  programs  for  this  analysis  were  written  in 
FOHTHAN  IV  for  execution  on  the  IBM  3033  attached  processor 
(System  370)  computer  installed  at  the  Naval  Postgraduate 
School.  The  programs  are  batch  oriented  since  the  volume  of 
data  available  prevented  storage  directly  on  the  user's  pri¬ 
vate  disk  space.  The  data  tapes  received  from  FMSO  were 
loaded  onto  the  system  mass  storage  device  and  when  required, 
all  or  portions  of  the  data  were  transferred  to  the  Virtual 
Machine  (VM)  facility  for  processing.  Two  programs  were 
written,  the  first  to  copy  the  data  to  make  it  accessible  to 
the  VM  system  and  the  second  to  conduct  the  actual  analysis. 

A  complete  listing  of  the  FORTRAN  source  code  is  available  in 
Appendix  F.  Each  program  or  subprogram  contains  a  documenta¬ 
tion  block  at  the  top  that  provides  a  complete  description  of 


the  program’s  purpose,  a  definition  of  variables  used  and  a 
list  of  user-written  subroutines  and  functions  required. 

The  random  number  generator  used  is  a  part  of  the 
LLRAHDOa  II  series  [Bef.  9]  developed  at  the  Naval  Postgradu 
ate  School.  The  programs  and  subprograms  were  thoroughly 
tested  both  independently  and  as  a  unit  using  known  data. 

The  results  were  verified  by  hand  held  calculator  and  in  all 
cases  agreed  with  the  program  output. 


V.  RESPLTS 


The  results  of  the  distributional  analysis  conducted  on 
the  ten  distribution  models  listed  in  Table  2  are  tabulated 
in  Appendix  B  for  1R  consumable  items.  Appendix  B  is  divided 
into  the  three  demand  classes  listed  in  Table  1,  and  then 
each  class  is  divided  by  distribution  and  percentile  within 
each  distribution.  The  results  for  the  2E  repairable  items 
are  presented  in  a  similar  format  in  Appendix  C. 

The  consistency  of  the  resampling  method  with  varying 
data  was  evaluated  by  testing  independent  sections  of  the 
consumable  and  repairable  data  sets  and  comparing  the 
results.  In  both  causes,  independent  subsamples  were  created 
by  sampling  every  third  item  but  changing  the  initial  item 
sampled  from  the  first  to  the  third  item.  The  analysis,  not 
presented  here,  showed  that  the  resampling  method  produced 
consistent  results  over  all  of  the  subsamples  tested.  The 
total  mean  square  error11  of  the  distributions  in  each  sub¬ 
sample  varied  by  no  more  than  20%  from  the  total  mean  square 
error  values  obtained  from  the  entire  sample  population  as 
listed  in  Appendix  B  and  Appendix  C.  The  relative  standing 


» ‘Total  mean  square  error  is  the  sum  of  the  mean  square 
errors  computed  at  the  75th,  80th,  85th,  90th  and  95th 
percentiles. 


of  each  distribution  was  not  altered  in  any  of  the 
sub sa spies. 

The  effect  of  the  nuab9r  of  pseudo  sample  repetitions  was 
also  studied.  The  results  proved  similar  while  varying  the 
number  of  repetitions  from  twenty  to  fifty.  The  standard 
deviation  of  the  percentile  estimates  decreased  as  the  number 
of  repetitions  increased.  The  total  mean  square  error  for 
each  trial  varied  randomly  but  remained  within  10H  of  the 
overall  values  for  each  distribution  with  more  repetitions 
producing  consistently  smaller  variations.  For  the  data 
analysis  production  runs,  forty  was  chosen  as  the  number  of 
resampling  repetitions  as  a  trade  off  between  computer  run 
time  and  minimizing  variance. 

A.  LOW-DEMAND  ITEMS 

A  summary  of  the  total  mean  squared  error  figures  for  the 
models  providing  the  best  fit  or  the  smallest  total  mean 
squared  error  is  provided  in  Table  3  for  the  low-demand 
items.  The  best  models  for  the  low  demand  items  all  had  very 
small  total  mean  squared  errors.  This  can  be  attributed  to 
the  very  low  probability  of  a  demand  ever  being  greater  than 
one  and  to  the  integer  sampling  plan  for  accumulating  the 
Bernoulli  trials.  For  the  integer  results  of  the  Poisson  and 
negative  binomial  distribut ions,  the  estimated  percentile 
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TABLE  3:  LOW-DEMAND  ITEMS, 

Total  Mean  Squared  Error  (10“4) 

Negative 

Poi  sson  binomial  Comegund 
Consumable  IB  5.  64  7.86  10.77 

Repairable  2B  5.  86  7.57  12.91 

will  closely  match  the  theoretical  percentile  unless  the  num¬ 
ber  of  Bernoulli  successes  is  less  than  the  expected  number 
of  successes  for  a  given  theoretical  percentile.  This  pro¬ 
duces  a  conservative  result  that  will  tend  to  overstate  the 
quantity  of  stock  required  by  one  unit  to  provide  a  stated 
level  of  protection,  but  in  the  case  of  the  small  quantities 
involved,  is  not  an  undesirable  preperty.  The  compound  dis¬ 
tributions,  as  a  group,  gave  the  same  total  mean  squared 
error  for  much  the  same  reasons.  The  probability  of  a  non¬ 
zero  demand  is  very  small  for  the  low-demand  items,  generally 
less  than  0.05,  thus  these  models  tend  to  predict  a  quantity 
of  zero  for  all  percentiles.  The  resampling  procedure  and 
the  sampling  plan  again  provide  a  conservative  result. 

B.  MEDIUM-DEMAND  ITEMS 

A  summary  of  the  total  mean  squared  error  figures  for  the 
models  providing  the  best  fit  or  the  smallest  total  mean 


squared  error  is  provided  in  Table  4  for  the  mediu a- demand 
items.  Here  again  the  predicted  demand  for  the  various  per¬ 
centiles  is  small.  The  resampling  technique  still  gives  the 

TABLE  4:  MEDIUM-DEMAND  ITEMS , 

Total  Bean  Squared  Error  (10~4> 

Negative  Bernoulli- 

biaoaial  Poisson  lognormal 

Consumable  IB  28.55  48.28  80.60 

Repairable  2B  29.14  44.  84  1  12.44 

edge  to  the  integer  distributions.  The  compound  distribu¬ 
tions  have  begun  to  show  some  differences  among  themselves, 
but  as  a  group  are  clearly  better  than  any  of  the  continuous 
simple  distributions.  The  magnitude  of  the  total  mean 
squared  error  is  significantly  larger  than  for  the  low-demand 
class  primarily  due  to  the  larger  allowable  total  demand 
range  defining  the  medium-demand  class. 

C.  HIGH-DEMAND  ITEMS 

The  summary  total  mean  squared  error  data  for  the  high- 
demand  items  are  listed  in  Table  5.  With  the  unlimited  range 
of  mean  demands  in  this  class,  the  higher  total  mean  squared 
errors  are  expected.  The  compound  models  have  come  into 
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TABLE  5:  HIGH-DEMAND  ITEMS, 

Total  Mean  Squared  Error  (10-*) 

Bernoulli-  Bernoulli- 

loqnor  na^.  exponential  Exponential 

Consuoable  IB  168.  25  188.90  298.21 

Repairable  2H  183.  77  221.80  289.01 

their  own  in  this  class  with  the  Bernoulli-lognormal  yielding 
a  total  mean  squared  error  less  than  one  half  that  of  the 
normal  distribution.  The  second  best  distribution,  the 
Bernoulli- exponential,  gives  a  total  mean  squared  error  only 
slightly  higher  than  the  Be rnoulli- lognormal. 


VI.  CONCLUSIONS  AND  RE  CO  MW  ENDATIONS 


A.  SUBBARY 

This  distributional  analysis,  conducted  using  a  cob- 
pletely  different  experimental  procedure  than  prior  analyses, 
provides  satisfyingly  similar  results  for  the  low  and  medium 
demand  categories  in  showing  that  the  Poisson  distribution 
gives  the  best  fit  for  the  low-demand  category  of  items  and 
that  the  negative  binomial  distribution  gives  the  best  fit 
for  the  medium-demand  category  of  items.  The  class  of  com¬ 
pound  distributions,  as  a  whole,  gave  good  fits  in  the  low- 
and  medium-demand  categories,  though  the  total  mean  squared 
error  was  two  to  four  times  as  large  as  that  of  the  Poisson 
or  negative  binomial  distribution. 

The  analysis  indicated  that  there  are  several  distribu¬ 
tions  which  provide  better  results  than  the  normal  distribu¬ 
tion  for  high-demand  items.  The  Bernoulli-lognormal 
distribution  consistently  provided  total  mean  squared  errors 
approximately  one  half  that  of  the  normal  distribution,  indi¬ 
cating  a  superior  fit  in  the  right-hand  tail  area  of  the  dis¬ 
tribution.  Other  distributions  giving  good  fits  were  the 
Bernoulli-exponential  and  the  standard  exponential,  both 
yielding  total  mean  squared  errors  less  than  that  of  the  nor¬ 
mal  distribution. 
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B.  RICOHHENDATIONS 


The  results  of  this  study  indicate  the  Navy  should  con¬ 
tinue  to  use  the  Poisson  and  negative  binomial  distributions 
as  the  models  for  the  low-  and  aediua-deaand  classes  respec¬ 
tively.  specifically,  the  Aviation  Supply  Office  should  re¬ 
evaluate  its  position  on  the  use  of  the  normal  distribution 
for  all  demand  classes  and  revert  to  the  Poisson  and  negative 
binoaial  distributions  for  low-  and  medium-demand  classes  as 
before.  The  normal  distribution  tends  to  inflate  the  stock 
required  for  low-  and  medium-demand  items  for  a  specified 
level  of  protection  resulting  in  increased  safety  levels  and 
excessive  dollar  investment. 

The  Navy  should  consider  replacing  the  normal  distribu¬ 
tion  model  used  for  high-demand  items  with  the  compound 
Bernoulli-lognormal  distribution.  If  the  Bernoulli-lognormal 
distribution  is  considered  computationally  too  difficult  for 
a  large  inventory  system,  the  Bernoulli-exponential  could  be 
used  in  its  place  with  little  loss  of  ef fectiveness. 


C.  ADDITIONAL  RESEARCH 

Pollow  on  research  in  the  following  areas  may  improve  and 
expand  upon  the  results  presented  above: 

1.  Optimize  the  divisions  of  demand  categories  (i.e.  low, 
medium,  high)  with  the  possible  use  of  the  probability 


of  a  non-zero  demand  as  an  element  in  the 
determination. 

Implement  the  Bernoulli-lognormal  model  in  FMSO's  5A 
simulator  to  establish  the  effect  of  the  change  on  the 
entire  supply  system. 

Investigate  the  implications  of  the  assumption  that  the 
demand  process  is  stationary. 

Evaluate  other  compound  distributions,  such  as  the 
Bernoulli -log-logistic,  which  is  particularly  appealing 
analytically. 


A££imx_A 


GLO??ABY  OF  T^RHS 

Program  Belated  Itea:  An  item  of  stock  whose  demand  can  be 
predicted  from  the  value  of  a  specific  Havy  program,  i.e., 
flying  hours,  steaming  hours,  etc. 

Beorder  Point:  The  on  hand  stock  quantity  that  when  reached, 
triggers  an  order  for  replenishment  of  stock  material.  The 
reorder  Point  is  the  expected  demand  during  procurement  lead- 
time  plus  the  safety  level. 

Bisk:  Probability  of  a  stockout  during  leadtime. 

Safety  Level:  The  quantity  of  material  which  is  required  to 
be  on  hand  to  permit  continued  operation  in  the  event  of 
minor  interuptions  on  normal  replenishment  or  unpredictable 
fluctuations  in  demand.  The  safety  level  determined  is 
structured  sc  as  to  minimize  time-weighted,  essentiality- 
weighted  requisitions  short. 


Stockout:  A  condition  that  exists  when  the  on  hand  inventory 

is  insufficient  to  fill  the  current  demand  requirements. 


System  Performance:  A  subjective  measure  maximized  when  th 
time-weighted,  essentiality-weighted  requisitions  short  is 
minimized.  Tiae  weighting  is  the  consideration  of  the  aver 
age  nuaber  of  days  delay  in  the  availability  of  material, 
essentiality-weight ing  is  the  consideration  of  the  relative 
essentiality  of  each  item.  Requisitions  short  are  requisi¬ 
tions  for  which  material  is  not  available. 


APPENDIX  B 

CONSUMABLE  ITEM  TABOLATED  RESULTS 
40  REPETITIONS  POR  EACH  OF  61  PSEUDO  SAMPLES 
LOW  DEMAND  RANGE:  0.0  <  D  <  1.0  PER  YEAR 


DISTRIBUTION;  POISSON 

THEORETICAL  PERCENTILES  (P) 


ESI  IS  HE  S: 

0.75 

0.85 

0*90 

MEAN 

0.  7500  00 

0.81  2500 

0.854006 

0.913924 

0 

STD  DEV 

0.008766 

0.007902 

0.007148 

0.005678 

0 

MSE 

0.0 

0.000156 

0.000018 

0.000290 

0 

TOTAL  MSE: 

5.64E- 04 

DISTRIBUTION: 

NEG ATI VE 

BINOMIAL 

THEORETICAL  PERCENTILES  (P) 

ESTIMATES: 

0.75 

0.80 

CK85 

0.90 

MEAN 

0.750000 

0.31  2500 

0.854006 

0.913924 

0 

STD  DEV 

0.008766 

0.007902 

0.007148 

0.005678 

0 

MSE 

0.0 

0.000156 

0.000018 

0.000290 

0 

0.95 

.956069 

.004149 

.000100 


0.95 
.948994 
. 004454 
.000322 


TOTAL  MSE 


7.86E-  04 


40  REPETITIONS  FOR  EACH  OF  61  PSEUDO  SAMPLES 
LOR  DEMAND  RANGE:  0.0  <  D  <  1.0  PER  YEAR 


DISTRIBUTION:  NORMAL 


MEAN 
STD  DEV 
HSE 

TOTAL  MSE: 


THEORETICAL  PERCENTILES  (P) 


0.75 

0*80 

0..85 

0..90 

0.95 

0.953221 

0.95  3221 

0.953221 

0.953221 

0.953221 

0.004275 

0.004275 

0.004275 

0.004275 

0.004275 

0.  042319 

0.024488 

0.01  1662 

0.003835 

0.001006 

833. 09E-04 


DISTRIBUTION:  EXPONENTIAL 

THEORETICAL  PERCENTILES  (P) 


1S2IHMIS: 

0.75 

0..85 

0._90 

0.95 

MEAN 

0.953221 

0.953221 

0.953221 

0.953221 

0.953221 

STD  DEV 

0.004275 

0.004275 

0.004275 

0.004275 

0.004275 

MSE 

0.042319 

0.024488 

0.011662 

0.003835 

0.001006 

TOTAL  HSE:  833.09E-04 


DISTRIBUTION:  LOGNORMAL 

THEORETICAL  PERCENTILES  (P) 


zsiizmz-- 

0*75 

0.80 

0. 85 

2i£0 

0.95 

MEAN 

0.  9784  04 

0.995246 

0.995246 

1.000000 

1.000000 

STD  DEV 

0.002943 

0.00  1392 

0.001392 

o 

. 

o 

0.0 

HSE 

0.052629 

0.03  8292 

0.021267 

0.009998 

0. 002500 

TOTAL  MSE:  1246.85E-04 
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40  REPETITIONS  FOR  E  ACH  OF  61  PSEUDO  SAMPLES 


LOW  DEMAND  RANGE:  0.0  <  D  <  1.0  PER  YEAR 


DISTRIBUTION:  LAPLACE 


£2XIMA2I£: 
MEAN 
STD  DEV 
MSE 

TOTAL  MSE: 


THEORETICAL  PERCENTILES  (P) 


fi*I5 

0  j_8  0 

0.85 

lo 

l» 

lu> 

lo 

0.95 

0.953221 

0.953221 

0.953221 

0.953221 

0.953221 

0.004275 

0.004275 

0.004275 

0.004275 

0.004275 

0.  042319 

0.024488 

0.011662 

0.003835 

0.001006 

833.09E-  04 


DISTRIBUTION:  LOGISTIC 

THEORETICAL  PERCENTILES  (P) 


ESTIMATES: 

0.75 

0.80 

0.85 

0.90 

0.95 

MEAN 

0.953221 

0.95  3221 

0.953221 

0.953221 

0.953221 

STD  DEV 

0.004275 

0.004275 

0.004275 

0.004275 

0.004275 

• 

MSE 

0.042319 

0.024488 

0.011662 

0.003835 

0.001006 

TOTAL  MSE:  833.09E-04 


DISTRIBUTION:  BERNOULLI-EX  PONENTI AL 

THEORETICAL  PERCENTILES  (P) 


0^75 

0^80 

0^85 

0.90 

0.95 

MEAN 

0.750000 

0.81  2500 

0.854006 

0.913924 

0.944939 

STD  DEV 

0.008766 

0.007902 

0.007148 

0.005678 

0.004618 

MSE 

o 

• 

o 

0.000156 

0.000018 

0.000290 

0.000613 

TOTAL  MSE: 

10.77E- 

04 
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40  REPETITIONS  FOR  EACH  OF  61  PSEUDO  SAMPLES 


snsii 

LOW 

DEMAND 

RANGE:  0.0  <  D  <  1 

.0  PER  YEAR 

DISTRIBUTION; 

BERNOU  LLI -LOGNORMAL 

THEORETICAL  PERCENTILES  (P) 

WH 

ESTIMATES: 

0*75 

0*80  0 t85 

2*20 

0.95 

.  * 

MEAN 

0.  750000 

0.31  2500  0.854006 

0.913924 

0.944939 

STD  DEV 

0.008766 

0.007902  0.007148 

0.005678 

0.004618 

i 

MSE 

0.0 

0.000156  0.0000 18 

0.000290 

0.000613 

- 

TOTAL  MSE: 

10.77E- 

04 

'  % 

DISTRIBUTION: 

BEENOU  LLI-LOGISTIC 

s 

>  ; 

» 

THEORETICAL  PERCENTILES  (P) 

. 

■i 

ESTIMATES: 

0*75 

0  .80  0.85 

2*20 

0.95 

■ 

MEAN 

0.7500  00 

0.81  2500  0.854006 

0.913924 

0.944939 

STD  DEV 

0. 008766 

0.007902  0.007148 

0.005678 

0.004618 

MSE 

0.0 

0.000156  0.000018 

0.000290 

0.000613 

TOTAL  MSE:  10.77E-04 


m 
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DISTRIBUTION:  POISSON 


THEORETICAL  PERCENTILES  (P) 


ESTIMATES: 

2*25 

2*22 

0  ..85 

0^90 

0.95 

MEAN 

0.  748131 

0.807719 

0.  843932 

0.893923 

0.925303 

STD  DEV 

0.002501 

0.00  2271 

0.002091 

0.001774 

0.001515 

MSE 

0.000138 

0.000398 

0.  000603 

0.001423 

0.002266 

TOTAL  USE:  48. 28  E-04 


DISTRIBUTION:  NEGATIVE  BINOtj^j, 

THEORETICAL  PERCENTILES  (P) 


222151112: 

2*21 

2*8 2 

2*25 

0.90 

0.95 

MEAN 

0.  746099 

0.805430 

0.843400 

0.900210 

0.943628 

STD  DEV 

0.002508 

0.00  2281 

0.002094 

0.001727 

0.001329 

MSE 

0.000332 

0.000522 

0.000589 

0.000807 

0.000606 

TOTAL  USE:  28.55E-04 


DISTRIBUTION:  NORMAL 

THEORETICAL  PERCENTILES  (P) 


ESTIMATES: 

2*25 

0  .80 

0.85 

0*22 

0.95 

MEAN 

0.  888030 

0.896699 

0.907882 

0.920152 

0.935573 

STD  DEV 

0.001817 

0.00 1754 

0.001666 

0.001562 

0.001415 

MSE 

0.024771 

0.01 4168 

0.007327 

0.003509 

0.002400 

TOTAL  MSE:  521.74E-04 
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40  REPETITIONS  FOR  EACH  OF  753  PSEUDO  SAMPLES 
MEDIUM  DEMAND  RANGE:  1.0  <  D  <  20.0  PER  I EAR 


DISTRIBUTION:  EXPONENTIAL 

THEORETICAL  PERCENTILES  (P) 


ESTIMATES: 

0*75 

0  .80 

0*85 

0.  90 

0*95 

MEAN 

0. 852981 

0.86 1003 

0.871591 

0.884311 

0.898352 

STD  DEV 

0.002040 

0.00  1993 

0.001928 

0.001843 

0.001741 

MSB 

0.020357 

0.01 1779 

0.006819 

0.004953 

0.006210 

TOTAL  MSE:  501.19E-04 


DISTRIBUTION:  LOGNORMAL 

THEORETICAL  PERCENTILES  (P) 


ESTIMATES: 

0.75 

2*80 

0.85 

0.90 

0.95 

MEAN 

0.953018 

0.96 2S52 

0.970638 

0.980623 

0.990731 

STD  DEV 

0.001219 

0.001094 

0.000973 

0.000794 

0.000552 

MSE 

0.043846 

0.02  8367 

0.015977 

0.007411 

0.002040 

TOTAL  MSE:  976.41E-04 


DISTRIBUTION:  LAPLACE 

THEORETICAL  PERCENTILES  (P) 


2*21 

0*82 

0*85 

0*90 

0.95 

MEAN 

0. 884254 

0.893767 

0.904733 

0.918110 

0.936078 

STD  DEV 

0.001843 

0.00  1775 

0.001692 

0.001580 

0.001409 

MSE 

0. 024079 

0.01 3865 

0.007191 

0.003573 

0.002363 

40  REPETITIONS  FOR  EACH  OF  753  PSEUDO  SAMPLES 


MEDIUM 

DEMAND 

RANGE:  1.0 

<  D  <  20.0 

PER  YEAR 

DISTRIBUTION: 

LOGISTIC 

THEORETICAL  PERCENTILES  (P) 

ESTIMATES: 

0*75 

0  .80 

0*85 

0*90 

0.95 

MEAN 

0.  8750  36 

0.886687 

0.896908 

0.913037 

0.934834 

STD  DEV 

0.001905 

0.00 1826 

0.001752 

0.001624 

0.001422 

MSE 

0.022871 

0.01 3193 

0.006955 

0.003690 

0.002447 

TOTAL  MSE: 

491. 56 E- 

04 

DISTRIBUTION: 

BERNOULLI -El  POHENTIAL 

THEORETICAL  PERCENTILES  (P) 

2*25 

0*80 

0*85 

0*90 

0.95 

MEAN 

0.738072 

0.795951 

0.840430 

0. 85i56  3 

0.950685 

STD  DEV 

0.002533 

0.00  2322 

0.  002110 

0.001755 

0.001248 

MSE 

0.002132 

0.00  2338 

0.002267 

0.002224 

0.001327 

TOTAL  MSE: 

102.38  E- 

04 

DISTRIBUTION:  BERNOULLI-LOGNORMAL 

THEORETICAL  PERCENTILES  (P) 


&XIfliX££: 

2*25 

2*80 

0*85 

2*22 

2*22 

MEAN 

0.  749344 

0.809159 

0.849391 

0.903886 

0.948263 

STD  DEV 

0.  002497 

0.002264 

0.002061 

0.001698 

0.001276 

MSE 

0.  001456 

0.00 1751 

0.001730 

0.001866 

0.001258 

TOTAL  MSE:  80.60E-04 
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40  REPETITIONS  FOR  EACH  OF  753  PSEUDO  SAMPLES 


MEDIUM  DEMAND  RANGE:  1.0  <  D  <  20.0  PER  YEAR 


DISTRIBUTION:  BERNOU  LLI-LO  GIST IC 

THEORETICAL  PERCENTILES  (P) 


ESTIMATES: 

0*75 

0^80 

0  ..85 

2*22 

0.95 

MEAN 

0.  7600  06 

0.82 1542 

0.864176 

0.916219 

0.955416 

STD  DEV 

0.002461 

0.00  2206 

0.001974 

0.001596 

0.001189 

MSB 

0.0022  56 

0.0024  99 

0.  0023  04 

0.002030 

0.001241 

TOTAL  MSE: 

103.30E- 

04 

1 

I 


40  REPETITIONS  FOR  EACH  OF  429  PSEUDO  SAMPLES 
HIGH  DEMAND  RANGE:  20. 0  <  D  PER  YEAR 


DISTRIBUTION:  POISSON 

THEORETICAL  PERCENTILES  (P) 


ESTIMATES: 

0.75 

2*80 

0.85 

2.90 

0.95 

MEAN 

0.689018 

0.73  0300 

0.753320 

0.780038 

0.800697 

STD  DEV 

0.003534 

0.00  3388 

0.003291 

0.003162 

0.003050 

MSE 

0.  024843 

0.027146 

0.033935 

0.042154 

0.051388 

TOTAL  MSE:  1794.66E-04 


DISTRIBUTION:  NEGATIVE  BINOMIAL 

THEORETICAL  PERCENTILES  (P) 


ESTIMATES: 

0.75 

2*80 

2*85 

0.90 

0,95 

MEAN 

0.724173 

0.747130 

0.791313 

0.848041 

0.896526 

STD  DEV 

0.003412 

0.00  3318 

0.003102 

0.002740 

0.002325 

MSE 

0.003047 

0.024647 

0.026816 

0.029773 

0.031637 

TOTAL  MSE:  1  159.20E-04 


40  REPETITIONS  FOR  EACH  OF  429  PSEUDO  SAMPLES 


HIGH  DEHAND  RANGE:  20.0  <  D  PER  TEAR 


DISTRIBUTION:  EIPONENTIAL 

THEORETICAL  PERCENTILES  (P) 


ESTIMATES: 

0*75 

0*8(2 

0.85 

2*90 

0.95 

MEAN 

0. 784610 

0.814275 

0.846685 

0.881207 

0.918702 

STD  DEV 

0.  003138 

0.002969 

0.002750 

0.002470 

0.002086 

MSB 

0.009534 

0.006885 

0.005219 

0.004384 

0.003800 

TOTAL  MSE: 

298.21E- 

04 

DISTRIBUTION:  LOGNORMAL 

THEORETICAL  PERCENTILES  (P) 


ESTIMATES: 

<2*25 

0  .80 

0*85 

0*90 

0.95 

MEAN 

0.855867 

0.886448 

0.918614 

0.949144 

0.978288 

STD  DEV 

0.002681 

0.00  2422 

0.002087 

0.001677 

0.001113 

MSE 

0.019028 

0.01  3031 

0.008113 

0.004389 

0.001639 

TOTAL  HSE:  462.00E-04 


DISTRIBUTION:  LAPLACE 

THEORETICAL  PERCENTILES  (P) 


ESTIMATES: 

0*25 

0 .80 

0*85 

0*90 

0.95 

MEAN 

0.  843416 

0.86  3600 

0.885612 

0.907473 

0.934956 

STD  DEV 

0.002774 

0.00  2620 

0.002430 

0.002212 

0.001883 

HSE 

0.015775 

0.00  9709 

0.005680 

0.003248 

0.002167 

TOTAL  HSE:  365.79E-04 
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40  REPETITIONS  POR  EACH  OP  429  PSEUDO  SAHPLES 
HIGH  DEHAND  RANGE:  20. 0  <  D  PER  TEAR 


DISTRIBUTION:  LOGISTIC 

THEORETICAL  PERCENTILES  (P) 


352231235: 

0x25 

0  .80 

0x85 

0.90 

0.95 

HEAN 

0.  822132 

0.846791 

0.871709 

0.9001  10 

0.933755 

STD  DEV 

0. 002919 

0.00  2750 

0.002553 

0.002289 

0.001899 

HSE 

0.  013721 

0.008841 

0.005633 

0.003523 

0.002259 

TOTAL  HSE:  339.76E-04 


DISTRIBUTION:  BERNOU LLI-EX PONENTIAL 

THEORETICAL  PERCENTILES  (P) 


352131135: 

0x25 

0x80 

0.85 

0x20 

0.95 

HEAN 

0.759754 

0.81  6600 

0.867492 

0.914225 

0.955101 

STD  DEV 

0.003261 

0.002954 

0.002588 

0.002138 

0.001581 

HSE 

0.  005828 

0.00  5034 

0.003977 

0.002618 

0.001432 

TOTAL  HSE:  188.90E-04 


DISTRIBUTION:  BERNOU  LLI  "LOGNORMAL 

THEORETICAL  PERCENTILES  (P) 


2x25 

0  .80 

0^8  5 

0x90 

0.95 

HEAN 

0.74  5078 

0.797614 

0.847684 

0.901963 

0.953670 

STD  DEV 

0.003327 

0.00  3067 

0.002743 

0.002270 

0. 001605 

HSE 

0.004768 

0.004403 

0.003632 

0.002613 

0.001409 

40  REPETITIONS  FOR  EACH  OF  429  PSEUDO  SAMPLES 
HIGH  DEMAND  RANGE:  20.0  <  D  PER  TEAR 


DISTRIBUTION:  BERNOULLI-LOGISTIC 

THEORETICAL  PERCENTILES  (P) 


ESTIMATES: 

2*25 

0.80 

0.85 

0..9Q 

0.95 

MEAN 

0.828539 

0.86  4722 

0.892571 

0.921326 

0  .947036 

STD  DEV 

0.002877 

0.00  2611 

0.002364 

0.002055 

0.001710 

MSE 

0.012517 

0.00  8862 

0.005485 

0.003020 

0.001538 

TOTAL  MSE 


314.23E-  04 


APPENDIX  C 


REPAIRABLE  IT  EM  y ABOLAT ED  RESULTS 
40  REPETITIONS  POR  EACH  OP  137  PSEODO  SAMPLES 
LON  DEMAND  RANGE:  0.0  <  D  <  1.0  PER  TEAR 


DISTRIBUTION:  gQjggON 


ESTIMATES: 
MEAN 
STD  DEV 
MSE 

TOTAL  MSE: 


THEORETICAL  PERCENTILES  (P) 


o*2S 

0.  750000 
0.005849 
0.0 

5.86E-  04 


0  .80 
0.81 2496 
0.005273 
0.000156 


0 185 
0.853613 
0.004775 
0.000019 


0.  90 

0.912594 

0.003815 

0.000311 


0*25 

0.955936 

0.002772 

0.000099 


DISTRIBUTION:  NEGATIVE  BINOMIAL 

THEORETICAL  PERCENTILES  (P) 


ESTIMATES: 

0.75 

0.80 

0. 85 

0&22 

0.95 

MEAN 

0.7500  00 

0.81 2496 

0.853613 

0.912594 

0.950753 

STD  DEV 

0.005849 

0.005273 

0.004775 

0.003815 

0.002923 

MSE 

o 

. 

o 

0.00  0156 

0.000019 

0.000311 

0.000270 

TOTAL  MSE 


7.57E-04 


h  - 

k 

40  REPETITIONS 

FOR  EACH  OF  137  PSEUDO  SABPLES 

m 

LOW 

DEBAND 

RANGE:  0.0 

<  D  <  1 

.0  PER  TEAR 

DISTRIBUTION 

:  NORMAL 

THEORETICAL  PERCENTILES  (P) 

S321&MMS: 

0.75 

0.80 

g.  85 

2*20 

0.95 

•> 

MEAN 

0.949164 

0.94  9164 

0.949164 

0.949 164 

0.949164 

- 

►  .  • 

STD  DEV 

0.002967 

0.002967 

0.002967 

0.002967 

0.002967 

m 

USE 

0.  040918 

0.023465 

0.01 1016 

0.003568 

0. 001118 

TOTAL  USE: 

800.  84 E- 

04 

vS 

DISTRIBUTION 

:  EXPONENTIAL 

\ 

»*..  •* 

• 

THEORETICAL  PERCENTILES  (P) 

i 

*  .  ■ 

ESTATES: 

0.75 

2*M 

0.85 

0.90 

0.95 

i 

* 

MEAN 

0.949164 

0.949164 

0.949164 

0.949164 

0.949164 

j 

» 

STD  DEV 

0.002967 

0 .00  2967 

0.002967 

0.002967 

0.002967 

'hag 

BSE 

0.  040918 

0.023465 

0.01  1016 

0.003568 

0.001118 

■ 

pm 

TOTAL  BSE: 

800.84E- 

04 

DISTRIBUTION 

:  iOS12M4i 

■4 

THEORETICAL  PERCENTILES  (P) 

EST  MATES: 

2*25 

2*12 

2*85 

2*22 

0.95 

.  . 

BEAN 

0.  978874 

0.99 1028 

0.991028 

1.000000 

1.000000 

STD  DEV 

0.001943 

0.001274 

0.001274 

0.0 

0.0 

: 

BSE 

0.052955 

0.036815 

0.020201 

0.009996 

0.002500 

4 

TOTAL  BSE: 

1  224.66  E- 

04 

, 

59 

4 

i 

u. •- 

- - »  -  *  -  ■■ 

J 

DISTRIBUTION:  LAPLACE 


THEORETICAL  PERCENTILES  (P) 


£2X£££IJ£: 

0,75 

0.80 

0 .85 

0.  90 

0.95 

MEAN 

0.949164 

0.949164 

0.949164 

0.949164 

0.949164 

STD  DEV 

0.002967 

0.00  2967 

0.002967 

0.002967 

0.002967 

HSE 

0.  040918 

0.023465 

0.011016 

0.003568 

0.001118 

TOTAL  MSE:  800.84E-04 


DISTRIBUTION:  LOGISTIC 

THEORETICAL  PERCENTILES  (P) 


ESTIMATES: 

0x15 

2x10 

2x85 

0.90 

0.95 

MEAN 

0.949164 

0.94  9164 

0.949164 

0.949164 

0.949164 

STD  DEV 

0. 002967 

0.002967 

0.002967 

0.002967 

0.002967 

HSE 

• 

0.  040918 

0.02  3465 

0.01 1016 

0.003568 

0.001118 

TOTAL  HSE:  800.84E-04 


DISTRIBUTION:  BERNOU LLI-EX PONENTIAL 

THEORETICAL  PERCENTILES  (P) 


ESTIMATES: 

2x25 

2x80 

2x85 

0x90 

0.95 

MEAN 

0.7500  00 

0.81 2496 

0.853613 

0.912594 

0.942635 

STD  DEV 

0.005849 

0.005273 

0.004775 

0.003815 

0.003141 

MSB 

0.0 

0.00  0156 

0.000019 

0.000311 

0.0 00804 

TOTAL  HSE:  12.91E-04 


60 


40  REPETITIONS  FOR  EACH  OF  137  PSEODO  SAMPLES 


LOW  DEMAND  RANGE:  0.0  <  D  <  1.0  PER  TEAS 


DISTRIBUTION:  BERNOD  LLI -LOGNORMAL 

THEORETICAL  PERCENTILES  (P) 


ESTIMATES: 

0.75 

ol 

2s.8£ 

0.90 

0.95 

MEAN 

0.7500  00 

0.81 2496 

0.853613 

0.912594 

0.942635 

STD  DEV 

0. 005849 

0.005273 

0.004775 

0.003815 

0.003141 

HSE 

TOTAL  MSE: 

0.0 

12. 91 E- 

0.000156 

04 

0.000019 

0.000311 

0.000804 

DISTRIBUTION:  BERNOU  LLI-LOGISTIC 

THEORETICAL  PERCENTILES  (P) 


ESTIMATES: 

2*25 

0.80 

2*25 

0 .90 

0.95 

MEAN 

0.7500  00 

0.81 2496 

0.853613 

0.912594 

0.942635 

STD  DEV 

0.005849 

0.00  5273 

0.004775 

0.003815 

0.003141 

MSE 

o 

i 

O 

0.000156 

0.000019 

0.000311 

0. 000804 

TOTAL  MSE: 

12.  91 B- 

04 

61 


40  REPETITIONS  POR  EACH  OP  1487  PSEUDO  SAMPLES 
MEDIUM  DEMAND  RANGE:  1.0  <  D  <  20.0  PER  I  EAR 


DISTRIBUTION:  POISSON 

THEORETICAL  PERCENTILES  (P) 


£5XIAAX2£: 

0*75 

0.80 

0.85 

2*22 

0*25 

MEAN 

0.747034 

0.806323 

0.843928 

0.894732 

0.929457 

STD  DEV 

0.001782 

0.00  1620 

0.001488 

0.001258 

0.001050 

HSE 

0.000238 

0.00  0507 

0.000582 

0.001321 

0.001836 

TOTAL  HSE:  44.84E-04 


DISTRIBUTION:  NEGATIVE  BINOM^L 

THEORETICAL  PERCENTILES  (P) 


ESTIMATES: 

o*I5 

2*22 

2*85 

0*20 

0.95 

MEAN 

0.744688 

0.804450 

0.843869 

0.900645 

0.944953 

STD  DEV 

0. 001788 

0.00 1626 

0.001488 

0.001227 

0.000935 

HSE 

0.000444 

0.00  0598 

0.000533 

0.000788 

0.000551 

TOTAL  HSE:  29.14E-04 


DISTRIBUTION:  NORMAL 

THEORETICAL  PERCENTILES  (P) 


ESTIMATES: 

2*12 

2*80 

2*85 

2*20 

0.95 

MEAN 

0.861474 

0.87  5334 

0.890316 

0.905783 

0.927155 

STD  DEV 

0.  001416 

0.00 1354 

0.001281 

0.001198 

0.001066 

MSB 

0. 019444 

0.01 1363 

0.006084 

0.003600 

0.002944 

TOTAL  HSE:  434.35E-04 


40  REPETITIONS  FOR  S  ACH  OF  1487  PSEUDO  SAMPLES 


MEDIUM  DEMAND  RANGE:  1,0  <  D  <  20.0  PER  YEAR 


DISTRIBUTION: 

£££oraiI!L 

THEORETICAL  PERCENTILES  (P) 

ESTIMATES: 

<2x25 

0,80 

2x85 

2x22 

0.95 

MEAN 

0. 823217 

0.834022 

0.845365 

0.864271 

0.886389 

STD  DEV 

0.001564 

0.00  1526 

0.001482 

0.001404 

0.001301 

MSE 

0.016435 

0.01 0140 

0.007542 

0.006691 

0.007993 

TOTAL  MSE: 

488.02E- 

04 

DISTRIBUTION:  LOGNORMAL 


I5XIfiA2S5: 

2x25 

THEORETICAL  PERCENTILES  (P) 

2x55  2x85  2x22 

0.95 

MEAN 

0.  946618 

0.955351 

0.964288 

0.975138 

0.987135 

STD  DEV 

0.000922 

0.000847 

0.000761 

0.000638 

0.000462 

MSE 

0.041761 

0.026454 

0.014803 

0.006692 

0.001890 

TOTAL  MSE: 

916. 00E- 

04 

DISTRIBUTION:  LAPLACE 

THEORETICAL  PERCENTILES  (P) 


ESTIMATES: 

2x25 

2x52 

2x55 

2x22 

0.95 

MEAN 

0.  856232 

0.870161 

0.885855 

0.902481 

0.928262 

STD  DEV 

0. 001439 

0.00  1378 

0.001304 

0.001216 

0.001058 

MSE 

0.  0188  23 

0.01 1026 

0.006113 

0.003736 

0.  002851 

TOTAL  MSE 


425.48E-04 


40  REPETITIONS  FOR  EACH  OF  1487  PSEUDO  SAMPLES 
MEDIUM  DEMAND  RANGE:  1.0  <  D  <  20.0  PER  YEAR 


DISTRIBUTION:  LOGISTIC 

THEORETICAL  PERCENTILES  (P) 


ESIIMAT^S: 

to 

L 

M 

0.80 

0.85 

0.90 

0.95 

MEAN 

0.  845409 

0.858922 

0.875792 

0.896303 

0.925823 

STD  DEV 

0.001482 

0.00 1427 

0.001352 

0.001250 

0.001075 

MSE 

0.0176  90 

0.01  06  03 

0.006269 

0.004086 

0.003093 

TOTAL  MSE:  417.42E-04 


DISTRIBUTION:  BERNOU LLI-EX PONENTIAL 

THEORETICAL  PERCENTILES  (P) 


ESTIMATES: 

0*75 

2^80 

0.85 

0*90 

0.95 

MEAN 

0.730065 

0.787917 

0.834508 

0.893026 

0.954267 

STD  DEV 

0.001820 

0.001676 

0.00 1524 

0.001267 

0.000857 

MSE 

0.003437 

0 . 00  36  08 

0.003449 

0.003021 

0.001465 

TOTAL  MSE:  149.79E-04 


DISTRIBUTION:  BERNOULLI -LOGNORMAL 

THEORETICAL  PERCENTILES  (P) 


ESTIMATES: 

2*25 

0*80 

2*85 

0*90 

0.95 

MEAN 

0.  750179 

0.808666 

0.852134 

0.904194 

0. 949477 

STD  DEV 

0.001775 

0.00  1613 

0.001455 

0.001207 

0.000898 

MSE 

0.  002226 

0.00  2681 

0.002600 

0.002327 

0.001410 

TOTAL  MSE:  112.44E-04 


40  REPETITION  FOR  EACH  OF  1487  PSEODO  SAMPLES 


MEDIUM  DEMAND  RANGE:  1.0  <  D  <  20.0  PER  YEAR 


DISTRIBUTION:  BBRNOfl LLI-LO  GIST IC 

THEORETICAL  PERCENTILES  (P) 


ESTIMATES: 

0x21 

2^80 

0 1 8  5 

0..90 

Q.  95 

MEAN 

0.7624  06 

0.82  2694 

0.866844 

0.915851 

0.953988 

STD  DE7 

0.001745 

0.00 1566 

0.001393 

0.001138 

0.000859 

MSB 

0.  003183 

0.00  3523 

0.  003107 

0.002468 

0.001379 

TOTAL  HSE 


136.  61 E-  04 


40  REPETITIONS  FOR  EACH  OF  356  PSEUDO  SAMPLES 


HIGH 

DEMAND  RANGE:  20.0 

<  D  PER 

YEAR 

DISTRIBUTION:  POISSON 

THEORETICAL  PERCENTILES  (P) 

ESTIMATES: 

2*75 

2*22 

0.85 

2*90 

0.95 

MEAN 

0.6  958  29 

0.735554 

0.762330 

0.791447 

0.821025 

STD  DEV 

0.003855 

0.00  36  96 

0.003567 

0.003405 

0.003212 

MSE 

0.022812 

0.026569 

0.032141 

0.038638 

0.044655 

TOTAL  MSE: 

1648. 16E-04 

DISTRIBUTION:  NEGATIVE  BINOMIAL 


2*15 

THEORETICAL  PERCENTILES  (P) 

0.80  0.85  0.90 

0.95 

MEAN 

0.722861 

0.755975 

0.799852 

0.857560 

0.904930 

STD  DEV 

0.003751 

0.00  3599 

0.003353 

0.002929 

0.002458 

MSE 

0.005285 

0.024724 

0.026795 

0.028699 

0.030113 

TOTAL  MSE: 

1156.16E- 

04 

DISTRIBUTION:  NORMAL 


551I3A US: 

2*25 

THEORETICAL  PERCENTILES  (P) 

2*22  2*25  0.90 

0.95 

MEAN 

0.  836334 

0.859810 

0.883941 

0.905503 

0.932035 

STD-  DEV 

0.003100 

0.002909 

0.002684 

0.002451 

0.  002109 

MSE 

0.014669 

0.009278 

0.005429 

0.003190 

0.002319 

TOTAL  MSE:  348.85E-04 


40  REPETITIONS  FOR  EACH  OF  356  PSEUDO  SAMPLES 


HIGH  DEHAND  RANGE:  20.0  <  D  PER  TEAR 


DISTRIBUTION:  EXPONENTIAL 

THEORETICAL  PERCENTILES  (P) 


ESTIMATES: 

Oils 

2i80 

0.85 

0.90 

0.95 

MEAN 

0.780940 

0.81 8187 

0.856974 

0.898153 

0.935255 

STD  DEV 

0. 003466 

0.003232 

0.002934 

0.002535 

0.002062 

MSB 

0.009578 

0.007296 

0.005229 

0.003830 

0.002968 

TOTAL  HSE:  289.01E-04 


DISTRIBUTION:  LOGNORMAL 

THEORETICAL  PERCENTILES  (P) 


ESTIMATES: 

0.75 

0 .80 

0..85 

0*22 

0.95 

MEAN 

0.819119 

0.85  5337 

0.895522 

0.933281 

0.971445 

STD  DEV 

0.003226 

0.00  2948 

0.002563 

0.002091 

0.001396 

MSE 

0.013003 

0.00  8704 

0.005985 

0.003420 

0.001454 

TOTAL  MSB:  325.65E-04 


DISTRIBUTION:  LAPLACE 

THEORETICAL  PERCENTILES  (P) 


ESTIMATES: 

2*11 

2  *80 

2*85 

0.90 

0.95 

MEAN 

0.  827652 

0.852088 

0.877334 

0.902360 

0.932870 

STD  DEV 

0.003165 

0.00  2975 

0.002749 

0.002487 

0.002097 

MSE 

0.  013659 

0.00  8798 

0.005346 

0.003302 

0.002267 

DISTRIBUTION:  LOGISTIC 


THEORETICAL  PERCENTILES  (P) 


1 SIUAUS: 

0*75 

0.80 

2*85 

0*2° 

0.95 

MEAN 

0.805812 

0.83  3134 

0.860882 

0.892308 

0.930748 

STD  DEV 

0.  003315 

0.003125 

0.002900 

0.002598 

0.002128 

MSE 

0.012232 

0.00  8347 

0.005666 

0.003803 

0.002397 

TOTAL  USE:  324.44B-04 


DISTRIBUTION:  BERNOU LLI-EX  PONENTIAL 

THEORETICAL  PERCENTILES  (P) 


ESIIHAIES: 

0.75 

0 .8Q 

2*85 

0*90 

0.95 

MEAN 

0.  770942 

0.829521 

0.877192 

0.921468 

0.960362 

STD  DEV 

0.003522 

0.003151 

0.002750 

0.002254 

0.001635 

HSE 

0.007108 

0.005916 

0.004548 

0.003092 

0.001516 

TOTAL  HSE:  221.80E-04 


DISTRIBUTION:  BERNOU  I.LI-LOGNORHAL 

THEORETICAL  PERCENTILES  (P) 


ESTIMATES: 

0*75 

0.80 

0*8  5 

0*90 

0.95 

MEAN 

0.  7461  08 

0.797150 

0.8471  97 

0.901039 

0.953617 

STD  DEV 

0.003647 

0.003370 

0.003015 

0.002502 

0.001762 

HSE 

0.005655 

0.004764 

0.003680 

0.002809 

0.001469 

TOTAL  HSE:  183.77E-04 
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40  REPETITIONS  FOR  EACH  OF  356  PSEUDO  SAMPLES 
HIGH  DEMAND  RANGE:  20.0  <  D  PER  TEAR 


DISTRIBUTION:  BERNOULLI-LOGISTIC 

THEORETICAL  PERCENTILES  (P) 


jSXiaiTjS: 

OjM 

0J35 

0^90 

0.95 

MEAN 

0.  829834 

0.86 1163 

0.886755 

0.914062 

0.941856 

STD  DEV 

0.003149 

0.00  2898 

0.002656 

0.002349 

0.001961 

MSB 

0.013074 

0.00  8963 

0.005495 

0.003143 

0.001722 

TOTAL  MSB 


323. 98E-04 


APPENDIX  D 


ESTIMATION  OF  LOGNORMAL  MEA S  AND  VARIANCE  FROM  UNTRANSFORMED 

m.h 

The  usaal  method  o£  estimating  the  sample  mean  and  vari¬ 
ance  of  a  random  sample  assumed  to  come  from  a  lognormal  dis¬ 
tribution  is  to  take  the  log  transform  of  the  data  and  then 
compute  the  standard  method  of  moments  estimates.  This  is 
time  consuming  on  a  computer  as  the  logarithm  must  be  taken 
for  each  observation  in  the  sample.  Formulae  to  compute  the 

mean  and  variance  of  the  un transformed  data  from  the  mean  and 

variance  of  the  transformed  data  can  be  derived  for  the  log¬ 
normal  distribution. 

If  Y  is  distributed  nor  mal(yy  ,o£)  ,  then  X  is  distributed 

ldgnornal  if  Y  *  ln(X).  The  mean  of  X  [Ref.  10]  is  expressed 


as: 


(20) 


and  the  variance  as: 


*  6 


Tfyy  +  TarVe °y 


-1  - 


(21) 


Equations  (20)  and  (21)  are  the  reverse  of  that  required  in 
this  case  since  and  a*  are  obtained  from  the  transformed 
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data.  What  is  needed  then,  is  to  solve  for  and  o*  in 
terms  of  jjx  and  oj.  Substituting  (20)  into  (21)  results  in: 


I*  a  ,.a 


which  when  solved  for  e  ,  taking  logs  and  simplifying,  gives 
the  solution  of  a*  as: 


/  2  2y 

y  \  , 


Substituting  (23)  back  into  the  log  of  (20)  ,  solving  for  uy 
and  simplifying  gives  the  solution  of  yy  as: 


^2  2 


Thus  the  mean  and > variance  cf  the  transformed  data  may  be 


estimated  from  the\  sample  mean  and  variance  of  the  untrans¬ 
formed  data  reducing  significantly  the  computational 


requirements. 
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C#*  XI MV  —  GIVEN  A  PROBABILITY ,  COHPUTES  THE  ABSCISSA  X  •** 

C**  VALUE  FROM  THE  DISTRIBUTION  INVERSE  FUNCTION 

C**  SAMPLE  —  GENERATES  PSEUDO  SAMPLES  AND  TESTS  AGAINST  ♦* 

THE  X  VALUE  COMPUTED  IN  XINV.  ACCUMULATES 
BERNOULLI  SUCCESS  AND  AVERAGE  MEAN  SQUARE  ERROR  ** 


r. 


CONDUCT  RESAMPLING  PROCEDURE  RETURNING  THE  NUMBER  OF 
BERNOULLI  SUCCESSES  AND  THE  AVERAGE  MEAN  SQUARED  ERROR  OF  THE 
MODEL. 


DO  60  J=1  ,  N  REPS 

CALL  SAMPLE  (IS  EED,  OBS, NBB,  I,  MODELS,  NPTILE  ,  X,  FX,  P  HAT,  MSE) 
CONTINUE 
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END 
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/' FILEDEF  ',  'TERMINAL ' ,  'DISK  ',  'OUTPUT 
•A1  *,  '(  ',  ' RECFM  ',  'FA 

•BLOCK  ',  *133  ',  'PERM  ',  'PILES 
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READ  (9  9,  110, EN  D=30)  IDSK  (I)  ,  DNAMEJ,  DTY PE$,  (OPT$(J)  ,J=1,5) 
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EXTRACT  CHARACTERS  ONE  AT  A  TINE  FROM  IVAL  AND  STORE  AS 
ALPHANUMERIC  CHARACTERS  IN  ARRAY  NUM 


IP  (ITEMP  .EQ.  0)  GOTO  20 
M  =  ITEMP 


CREATE  OUTPUT  ARRAY 


STBS  =  B$ 

IP  (IPRT  .GE.  30)  WRITE  (ICNOOT,  220)  IVAL,STR $ 

220  FORMAT ('  HOMERIC  TALO  E  =  ',18,*  ALPHANUMERIC  VALUE 
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IPBT,  ICONIN,  ICNOUT,  IPRINT 
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SUBROUTINE  OUTPUT  (MODEL, NPT ILE, NREPS,FX, IC,PHAT,MSE,CAT) 
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REPETITIONS  FOR  EACH  PERCENTILE  AND 
MODEL  TESTED  *» 

TMSE  —  TOTAL  MEAN  SQUARE  ERROR  SUMMED  OVER  PERCENTILES  $3 

FOR  EACH  DEMAND  CATEGORY  AND  MODEL  0$ 

TRIALS  —  TOTAL  NUMBER  QF  REPETITIONS.  USED  TO  NORMALIZE 


C**  PINAL  FIGURES.  COMPUTED  AS  COUNT  OF  LIN  I  ITEMS 

C«#  PER  DEMAND  CATEGORY  TIMES  NREPS 
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